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We report the long-time nonlinear dynamical evolution of ultracold atomic gases in the P-band 
of an optical lattice. A Bose-Einstein condensate (BEG) is fast and efficiently loaded into the P- 
band at zero quasi-momentum with a non-adiabatic shortcut method. Eor the first one and half 
milliseconds, these momentum states undergo oscillations due to coherent superposition of different 
bands, which are followed by oscillations up to 60ms of a much longer period. Our analysis shows the 
dephasing from the nonlinear interaction is very conducive to the long-period oscillations induced 
by the variable force due to the harmonic confinement. 

PACS numbers: 67.85.-d;03.75.Lm;03.75.Hh; 37.10.Jk 


I. INTRODUCTION 

The occupation of high orbital quantum state in solid- 
state system is known to play a crucial role to the un¬ 
derstanding of a lot of strongly correlated systems, such 
as high-temperature superconductivity and semiconduc¬ 
tor [l|. Ultracold atoms in an optical lattice have been 
used to simulate quantum many-body systems in the 
ground state Q. Recently, there has been increasing ef¬ 
forts to study ultracold atomic gases in optical lattices 
at high orbitals Q . This kind of study can not only help 
to understand high orbital quantum state in solid-state 
system but also provide a chance to study new physics 
which is completely beyond that of the solid-state sys¬ 
tem, such as the emergence of the quantum coherence in 
the first excited Bloch band of an optical lattice Q . High 
orbital physics in cold atomic gas has attracted a lot of 
attention; there are many theoretical predictions, such as 
supersolids [5| and other novel phases [6|-[8|. 

The long time dynamics of a Bose-Einstein condensate 
(BEG) in the quantum state of the S-band has been in¬ 
tensively studied and many interesting phenomena, such 
as coherent oscillations [3, [ifl nonlinear self-trapping 
pm were observed. However, it is much more chal¬ 
lenging to observe the long-time dynamics of a BEG in 
high orbitals because of the difficulty to prepare the cold 
atoms in such states and maintain their quantum coher¬ 
ence. For the condensate with a wide quasi-momentum 
distribution in P-band Q , the quantum dynamics for the 
holding time being shorter than 1.3 ms were studied. 

Here we experimentally prepare a BEG in the P- 
band of an optical lattice using a non-adiabatic short¬ 
cut method and study its long time quantum dynamics. 
The initial BEG in the P-band with finite size at quasi- 
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momentum q = 0 has a local maximum in energy so that 
this state is unstable when the harmonic confinement and 
interatomic interaction are considered. For the first one 
and a half milliseconds, the momentum states undergo 
fast oscillations. After a brief intermediate interval, the 
BEG starts a different type of oscillations which have 
a much longer period 14.9ms) and last up to 60ms. 
We find that these long-period oscillations are due to the 
presence of an additional harmonic trap, and the dephas¬ 
ing induced by the nonlinear interaction is greatly helpful 
for this oscillation behavior. This counterintuitive fea¬ 
tures and long time quantum coherence are possible to 
study the quantum thermalization. 

II. PREPARATION OF P-BAND QUANTUM 
STATE WITH A SHORTCUT METHOD 

In our experiment, we first prepare a pure BEG of 
about 2 X 10^ ^^Rb atoms in a hybrid trap which is 
formed by overlapping a single-beam optical dipole trap 
with wave length 1064nm and a quadrapole magnetic 
trap. The resulting potential has harmonic trapping fre¬ 
quencies {uJx,ojy,uJz) = 27 r X (28, 55, 65)Hz, respectively. 
Our lattice is produced by a standing wave created by 
two counter-propagating laser beams with lattice con¬ 
stant a = A/2 = 426 nm along the x-direction. 

We then load the BEG into the P-band by applying a 
series of pulsed optical lattices within tens of microsec¬ 
onds |13l-[lq. There are many other methods prepar¬ 
ing a BEG in the P-band quantum state Blilil, but 
normally the loading time is tens of milliseconds. Our 
method has the merit of rapid generation of the quan¬ 
tum state in the P-band. More importantly, this method 
allows us to generate the desired P-band with very nar¬ 
row quasi-momentum width around q = 0, which is cru¬ 
cial to observe the long-time quantum dynamics. Most 
recently, this method is proposed to study anisotropic 
two-dimensional Bose gas . 
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FIG. 1: (al-a5) The first row: the measured momentum dis¬ 
tributions of the BEG in the P-band at holding different times 
(0ms, 2ms, 5ms, 7ms, 30ms), and the white rectangles are re¬ 
gion for us to calculate proportion of different momentum 
states; the second row: schematic illustration of the corre¬ 
sponding population distributions in the Bloch band, (b) 
Population oscillations around momenta Ohk (blue stars) and 
—2hk (red dots) with t < 1.5ms. Solid lines are from numeri¬ 
cal simulations. The shade is to highlight the slow oscillating 
envelop of the rapid oscillations, (c) Population oscillations 
around momenta Ohk (blue stars), —hk (black diamonds), and 
—2hk (red dots) with t > 2ms. The oscillation period is about 
14.9ms with Vo = 5Er. 


We use two acousto-optic modulators to form our 
designed pulse sequences with the frequency difference 
5cjJ = 182.5MHz which corresponds to a phase shift be¬ 
tween two pulses series by 37r/4. This phase shift is cru¬ 
cial to change the parity of the ground state and thus 
prepare the state in the P-band. Four special pulses with 
a chosen depth Vq are used to prepare the state in the 
P-band around q = 0. The details of the pulse series are 
shown in the Supplementary Material. 

For the quantum state in the S-band with q = 0^ the 
momentum distribution has peaks around Ohk and ±2fi/c 
with k = 27 t/X [mil. For the quantum sate in the 
P-band with q = 0^ because its wave function has the 
odd parity 2 p{—x) = —'^(x), the momentum distribution 
equals to zero at Ohk and has significant peaks at ±2fi/c, 
as shown in Fig. l(al), where two dominating peaks at 
±2hk are clearly seen. 


III. DYNAMICAL EVOLUTION OF THE 
QUANTUM STATE IN P-BAND 


To study the dynamical evolution after a BEG is pre¬ 
pared in the quantum state of the P-band, we hold the 
condensate for time t before measuring the momentum 


distribution of the ultracold atoms after a time of flight 
(TOF) of 28ms. The momentum distributions of the 
BEG at t =0ms, 2ms, 5ms, 7ms, 30ms are shown in 
the first row of Fig. l(al, a2, a3, a4, a5). An evolu¬ 
tion is clearly observed. The corresponding evolution in 
the Bloch bands is schematically illustrated in the second 
row. The depth of the optical lattice is Vb = with 
Er being the atomic recoil energy. 

The momentum distributions are analyzed by comput¬ 
ing the normalized atom populations around momentum 
states \lhk) {I = 0, —1, —2) as Wi{t) = \{lhk\'ip{t))\‘^ = 
Ni{t)/N. Here N is the total atom number and Ni{t) is 
the atomic number around \lhk). In the TOF images we 
choose 95.2//m as the width to determine the atom num¬ 
ber Ni. This is the width of a condensate without any 
quantum manipulation and after the free expansion of the 
same TOF. We find W’s oscillation with time as shown 
in Fig. l(b, c). In Fig. 1(b), there are rapid oscillations 
in W- 2 {t) and Wo{t). These rapid oscillations disappear 
around t = 1.5ms. After a short transition time, a differ¬ 
ent type of oscillations begin to emerge around t =2ms. 
These oscillations can be observed up to 60ms and have a 
much longer period of 14.9ms. We will show how the in¬ 
teraction would affect this long-period oscillation, which 
was not observed in the previous work Q. 


IV. SHORT-PERIOD OSGILLATIONS 

The short-period oscillations shown in Fig. 1(b) are 
beating signal due to coherent superposition of different 
bands. In real experiment, it is impossible to prepare a 
quantum state completely in the P-band. What we pre¬ 
pared is in fact a superposition of s, p and d bands with 
negligible higher band. Due to difference in timescale of 
two oscillations and uncertainty in experiment, the ex¬ 
perimental data is shown in the figure comparing with 
the peripheral contour of the beating signal from the¬ 
oretical simulation with the choice of the initial quan¬ 
tum state |'0(t = 0)) = VK9\p,q = 0) + V0.05|d, q = 
0) + A/aM|s,g = 0). 

With this understanding, the periods of the oscilla¬ 
tions in Fig. 1(b) are determined by the band gaps. As 
there are three band gaps between the s, p, and d bands, 
there should be three oscillating periods, Tgp, Tgd, and 
Tpd- Around Ohk^ the atoms are mostly from the s and d 
bands due to the odd parity of the P-band. As a result, 
we expect only one oscillation period Tgd = 60.0/is. This 
is indeed what we see for Wo in Fig. 1(b). Around —2hk^ 
all three bands have significant contributions. Because 
the populations of the s and d bands are both small, the 
contribution to W -2 related to period Tgd is of second- 
order and negligible. As a result, only two periods are 
expected for W_ 2 , Tgp and Tpd- The oscillations of W -2 
in Fig. 1(b) clearly have two different periods, rapid os¬ 
cillations enveloped by slow oscillations. This is consis¬ 
tent with the fact that Tpd is much larger than Tgp. For 
Vb = SE’r we have Tgp = 69.5/is and Tpd = 465.5/is. 
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The short-period oscillations in Fig. 1(b) have a slow 
decay. This decay is mainly due to the decay of popula¬ 
tion on P-band induced by collision between atoms 
which is already considered in theoretical results in 
Fig. 1(b) (solid lines). For our system with atomic num¬ 
ber 2 X 10^ and density of 6.4 x 10^^cm“^ for Vb = SF’r, 
the calculated decay time is 1.3ms. It is coincident to the 
fitting result (1ms) based on the experimental data, and 
the decay rate in experiment is larger which can be at¬ 
tributed to other mechanism including dephasing of the 
atom cloud. Considering its finite size, the condensate 
has a narrow quasi-momentum distribution around q = 0 
with a width of O.Odhk. Our numerical simulation finds 
that the calculated decay time would increase by about 
1.9%, and numerical simulation shows the influence of 
this distribution to oscillation frequency is also negligi¬ 
ble. 


V. LONG-PERIOD OSCILLATIONS IN THE 
P-BAND 

The rapid oscillations disappear at about 1.5ms. After 
a featureless short dynamical evolution, a different type 
of oscillations begin to emerge at around 2ms. As shown 
in Fig. 1(c), these oscillations have a much longer period 
of about 15ms and their oscillating amplitude subsides 
gradually and becomes zero around 60ms after five cycles. 
This is different to the oscillation period of the harmonic 
trap 35ms. If we consider this evolution based on the 
semi-classical model with an effective potential, where 
the atoms exerted a variable force due to the harmonic 
trap and quasi-momentum dependent effective mass, as 
shown in the Supplementary Material, we see that the 
period is identical to the experimental result. Different 
to the Bloch oscillation for atoms on S-band [13, [H, llj ? 
however, our analysis shows that these long-period oscil¬ 
lations for atoms on P-band are greatly affected by the 
nonlinear dephasing, and this semi-classical model can 
not fully coincident to the experimental observation. 

We now look at the experiments in detail. Shown 
in Fig. [2fa) are the momentum distributions measured 
experimentally for holding times up to 15ms, which is 
roughly the first cycle of the oscillations. At t = 0ms, 
the momentum mainly distributes around ±2fi/c with 
q = 0, which is the signature of the initial state in the 
P-band. At t = 3ms, the absolute value of the peak 
quasi-momentum becomes smaller than 2hk. This is due 
to that the kinetic energy is transferred into the poten¬ 
tial energy by the harmonic trap. At t = 5ms, the atoms 
spread to the full regime between —2hk and -\-2hk. At 
t = 7ms, the atoms concentrate around Tfi/c, an indi¬ 
cation of the quantum state has moved to the Brillouin 
zone edge of the P-band (Fig. I(a4)). Around 15ms, most 
of the atoms have moved back to positions around ±2hk. 
Overall, an oscillating pattern is clearly observed in these 
momentum distributions. 

The oscillation pattern can be well captured by our 
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FIG. 2: The two-dimensional momentum distributions of the 
first cycle are shown in (a) and (b) for experiment and the¬ 
oretical simulation, respectively. After an integration along 
the vertical direction, (c) gives the one-dimensional distribu¬ 
tion for experiment (red dotted line) and numerical simulation 
(blue dashed line). The black line is obtained from the semi- 
classical model without considering interatomic interaction. 


simulation with the Gross-Pitaevskii equation (GPE) as 
seen in Fig. [2](b), where the simulation is to begin with 
the ground state in harmonic trap, and carry out with the 
same experimental process. At the center of the momen¬ 
tum distribution, the experiment shows significant popu¬ 
lation, which is due to the presence of incoherent atoms. 
To show it more clearly, we have integrated out the verti¬ 
cal direction and compared the experimental and numer¬ 
ical results in Fig. [2fc). We see that the peak positions 
are seen clearly matching very well. The experimental 
results (red dotted line) agree with the theoretical sim¬ 
ulation (blue dashed line), where the black line is based 
on the semi-classical model. 


VI. EFFECTS OF THE RANDOM PHASE 
BETWEEN NEIGHBORING LATTICE SITES ON 
THE LONG-PERIOD OSCILLATION 

The observed time evolution (black dots with error 
bar) of the proportion with momentum Hk is given in 
Fig. [3](a) for Vb = 5E^, where we can distinguish five 
cycles of oscillation. In the figure we have extracted in¬ 
coherent atoms at the center, and correspondingly the 
numerical simulation is multiplied by a decay factor. For 
bEr, the theoretical simulation with interatomic inter¬ 
action (red solid curve) agrees perfectly with the exper- 
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imental result. It is quite surprising that the simula¬ 
tion without interatomic interaction (blue dashed curve) 
doesn’t show clear periodic behavior, which differs sig¬ 
nificantly with the experiment. Both in experiment and 
theoretical simulation, the interaction has an effect of 
stabilizing the oscillation. This is counterintuitive by 
noticing that (i) the repulsive interaction has the effect 
of widening the wave packet and thus it seems that it 
should destroy the oscillation; (ii) the interaction usually 
plays a role of destroying the quantum coherence because 
of the nonlinear interaction. 

To solve this puzzle, we consider the following evolu¬ 
tion of the order parameter of the condensate. 
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Here Nj is the number of atoms initially in the jth lat¬ 
tice site, and '0j(x, t) = ajn{t)Wn{x) with Wn{x) the 
Wannier wave function located in the nth lattice site and 
OLjn{t) the corresponding coefficient. For atoms in jth 
site, initial situation is o:jn{^) = During evolution 
of the single site wave function, population would trans¬ 
fer from the jth site to different site n, and ajn{t) can be 
rewritten as Despite a smoothly varying 

of the external harmonic potential, the contribution of 
both the kinetic energy and potential energy to the phase 
(j)jn{t) is the same for every lattice site. However, the 
phase due to the interatomic interaction is proportional 
to each site’s atom number which could be quite different 
for neighboring lattice sites. As shown in Fig. [3l^c), the 
phase difference between neighboring lattice sites appears 
completely random because of interatomic interaction at 
t = 5ms. In a sense, the relative phase between neigh¬ 
boring sites is pseudorandom because it appears random 
only after the relative phase is confined between — tt and 


In this case, when the density distribution or momen¬ 
tum distribution is considered, the interference term be¬ 
tween different ajn{t) could be omitted, i.e. 
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y ] 0!jn{t)Wn{x) 
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With the above approximate expression, we solve every 
'ipj{x,t) by the discrete nonlinear Schrodinger equation 
(DNLSE) 0, and then get the momentum distribution, 
as shown with the green dotted lines in Figs. lH^a) and 
[3{b). In this situation, our numerical calculations show 
periodic behavior with each peak shift ahead a little com¬ 
paring with experiment, however, the period of the time 
agrees very well with the experimental result. In the 
Supplementary Material, we give the evolution of the 
|'0j(x, which shows that there is long time period be¬ 
havior in For the condensate in the P-band, 



FIG. 3: The population with hk versus time t for bEr (a) 
and lOEr (b). The experimental results, theoretical simu¬ 
lations with and without the interaction correspond to the 
black dots with error bar, red solid curves and blue dashed 
curves, respectively, and the green dotted curves are the re¬ 
sult of DNLSE. (c) The phase differences between neighboring 
lattice site j for Vo = bEr and t = 5ms. (d) The dotted points 
are the measured decay time of the oscillation amplitude, and 
the cross points are for the dynamical instability ratio. 


because in every it has the peaks around —2hk 

and 2hk^ the wave packet in every lattice site will split 
into two parts, propagate in different directions and os¬ 
cillate due to the presence of the harmonic trap. It is nat¬ 
ural that the incoherent superposition in Eq. (|2]) will also 
have the long time periodic oscillation. For Vq = bEr the 
GPE agrees better with the experiment than the method 
of DNLSE and Eq.([2|), while here the DNLSE simulation 
shows the physical picture more clearly. 

To verify further the role of random phase in the ap¬ 
pearance of the long-time period behavior, we add this 
random phase artificially for noninteracting case and 
solve numerically the Schrodinger equation. We find that 
there is perfect periodic behavior in IFi(t), although the 
period is a little different from the interaction case. 

In Fig. [3Kb), we give the experimental result for Vq = 
10Er. We see that the amplitude of the oscillation for 
10Er is much smaller than bEr. For lOE'r, there is a much 
stronger damping which is due to the highly suppressed 
tunneling between neighboring sites. In addition, there 
is no clear periodic behavior in this case. The simulation 
with GPE doesn’t agree well with the experimental re¬ 
sult. This may due to the fact that at 10E^, the particle 
number is highly squeezed so that there is a breakdown of 
the mean-field theory, which induces the disappearance 
of the interference between It seems that the 

theoretical simulation with GPE for bE^ agrees better 
with the experiment than the method with DNLSE. 
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VII. DISCUSSION AND SUMMARY 

In the long time evolution, the oscillation amplitude 
and damping are greatly affected by the dynamical insta¬ 
bility due to the interaction between atoms. We calcu¬ 
lated the dynamical instability [mil of P-band at lat¬ 
tice depth varying form 5 to lOE’r- In Fig.[3fd) the cross 
points show the calculated dynamical instability with the 
ratio of different lattice depth and bEr. The solid circles 
give the lifetime of oscillation by fitting the experimental 
data with exponential decay. As expected, the increasing 
of Vq will lead to larger dynamical instability ratio, and 
lead to the suppression of the oscillation. 

In summary, we demonstrate an effective way of load¬ 
ing atoms directly into P-band using two group of lat¬ 
tice pulses with different phases to break the parity of 
the Hamiltonian. Our preparation of the P-band with 
g = 0 is different from Ref. Q where the condensate 
has a distribution in the regime between q = —hk and 
hk. In our work, the evolution of the condensate in the 
P-band is studied experimentally. The short-period os¬ 
cillation due to coherent superposition of different bands 
and long-period oscillation reflecting the random rela¬ 
tive phase between neighboring lattice sites are observed 
simultaneously. The present experiment paves the way 
to study the long time dynamical evolution of the high 
orbital physics for other novel quantum state such as fer- 
monic superffuid, molecular condensate and condensate 
with special configuration of optical lattice. 
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Appendix A: Loading method 

The Hamiltonian of a single-atom in a one-dimensional 
optical lattice can be expressed as: 

Hx = P^/2m -h Vb cos^ {kx ), (Al) 

where Vb is the depth of the optical lattice, which is in 
unit of Er^ Er = {hk)‘^l2m is the recoil energy of the 
atom, px is the momentum of the atom along x-direction, 
m is the atom mass, /c = 7r/a is the wave vector and a is 
the lattice constant. The eigenstate of the Hamiltonian 
Eq. (|Aip is the Bloch state |n,g). 

Since the system is periodic, and our loading method 
conserves the quasi-momentum, we can choose the plane 
waves \2lhk^q) with quasi-momentum g as a set of com¬ 
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FIG. 4: The shortcut loading process of the atoms into P-band 
in Id optical lattice, (al) The schematic diagram of atoms in 
P-band. (a2) Loading time sequences with a different phase 
between the first series of pulses from 0 to and the second 
from ti to t 2 . t is the holding time in the optical lattice, {bl- 
bS) illustrate the superposition coefficient ci right before ti 
and c'l at the moment right after ti and at t 2 , respectively. 
(cl-c3) illustrate the corresponding population distribution in 
the Bloch band. 


plete bases, i.e., 

\ip) = ^ci\2lhk + q), (A2) 

I 

where q is the superposition coefficient, I = 0,±1,±2... 
is an index number, g G [—/c, k) is the quasi-momentum, 
and \2lhk-\-q) = is the eigenstate ofpx- 

On the other hand, the state lip) can also be expressed 
in the Bloch state |n,g) as \ 2 p) = q\ 2 p)\n, q)^ where 

n = s^p^d... is the band index. The aim state {ipa) is the 
first excited state at g = 0 , as shown in Fig. lU^al). 

Our shortcut method includes a series of designed 
standing wave pulses. The atom state \ 2 p) = 

n ■ Up.(tp^) ■ \ipo), where Up^{tp.) and 

j = l 

are the time evolution operator and free evolution oper¬ 
ator for the jth pulse with duration tp. and tf.^ respec¬ 
tively. The atom experiences spatial period potential in 
the pulse duration, and the different momentum states 
can gain different phase factors between the interval. 

Our preparation process consists of two series of pulses 
as shown in Fig. Ill^a2), where we need optimize the fi¬ 
delity |('0a|'0)P between the aim state \2pa) and \ 2 p) to its 
maximum by changing the parameters of pulse sequence. 
In the first series of pulses from 0 to G, the atom expe¬ 
riences spatial potential Vb cos^ {xn/a). For the second 
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series pulses from ti to t 2 the atom experiences poten¬ 
tial Vb cos^ (x7r/a-f 37r/4). The coefficients q (cJ) and 
the distribution in the Bloch band are shown in Fig.lU^b) 
and (c), respectively. 

For a pure Bloch state or a superposition state with 
different Bloch state at = 0, the parity can be given as 
P = "^\ci — c_/p/4, where P = 1 stands for a state with 
i 

odd parity and P = 0 even. 

At time ti, all the components in the parity P would 
satisfy q — C-i = 0, as shown in Fig. 111^61), and the 
atom is distributed in the s,d,g... bands, as in Fig.lU^cl). 
However, from the view of the second series of pulses, by 
the lattice shift, the coefficient q would be appended with 
a phase according to /, as c[|/) = |/), and the 

relation between the coefficients become — (—= 
0. In our loading process, the first series of pulses ensure 
that coefficient cj with even I is zero, and thus the parity 
of state can be completely changed as shown in Fig. 01^62). 
Its energy band distribution is in the p,f,.., as drawn in 
Fig. |4l^c2). At time ^ 2 , as shown in Figs. 111^63) and (c3), 
with another two pulses conserving the parity, only P- 
band state is populated. 

Appendix B: Effective potential model for p band 
dynamics 

Here we consider a semi-classical model to understand 
the periodic dynamics in the P-band. In the optical 
lattice, the effective mass is m*{q) = h{dv/dq)~^ = 
h‘^{d‘^Ep{q)/dq‘^)~^, where V is the group velocity and 
Ep{q) is the energy of the P-state for different q. In 
Fig. [5fa), the ratio between m and m* is shown for 
Vo = 5Pr and lOP^, respectively. We see that the effec¬ 
tive mass changes sign at the q = ±.^.21hk and ±0.36hk 
for Vo = 5Pr and lOPr, respectively. 

With the semi-classical model, we simplify the problem 
by considering a quasi-particle with q = 0 and study the 
dynamical process in the harmonic trap. When the ef¬ 
fective mass is considered, the acceleration is given as 
a{x) = [ 23 , where E{x) = —d{^muj‘^x‘^)fdx = 

—muj‘^x is the force driven by the harmonic trap with¬ 
out the optical lattice. Using ^ = 
and the initial condition, we get q, m*, x and E at time t. 
Then by expressing all the parameters as function of x, 
the dynamical evolution can be described by an effective 
potential: 

In Fig.[5{b), we give the numerical results of the effective 
potential. This means that the external confinement acts 
as a repulsive potential [13, [HI and the atoms at trap 
center would be divided into two symmetrical parts for 
the different momentum states. The minimums of the 
potential correspond to the points where the sign of the 
effective mass changes from the negative to positive. The 



q(7ik) 



FIG. 5: The effective mass for atoms in P-band (a) and its 
effective potential in the presence of a harmonic trap (b) for 
5Pr(red solid line) and 10Pr(blue dashed line). The dotted 
lines correspond to the points where the sign of the effective 
mass changes. 


initial atoms in the P-band at g' = 0 and x = ±5/rm near 
the center of the effective potential, would move from the 
center to the edge (as points denoted by time 1, 3, 5, 7ms) 
and then come back (as points by time 9,11,13,15ms) 
and form an oscillation. 


Appendix C: Simulation with single site DNLSE 

Here we demonstrate the long-period oscillation from 
the random phase between neighboring lattice sites in de¬ 
tail. As shown before, the simulation without interatomic 
interaction doesn’t have a perfect oscillation. However, 
when the interference terms between different sites are 
omitted, we can simulate the evolution of each single site 
based on the DNLSE Q, and the population rij would 
evolute with time as shown in Figs. [6] and 0 for bEr and 
lOEr lattice depth, respectively. For SE'r, atoms in differ¬ 
ent lattice sites would be distributed more concentrated 
during the oscillation, while for 10Er the atoms would 
disperse more rapidly into different lattice sites. 

According to the evolution of the wave function in the 
jth lattice site t) of the optical lattice and the har¬ 
monic trap, for each site’s evolution, there is a perfect 
oscillation for the ratio of atoms with momentum hk ver¬ 
sus the holding time in the potential. The initial state 
T(x,t = 0) is the superposition of Wannier functions in 
different lattice sites. Considering the fact that the ran- 
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domization between neighboring lattice sites arises from 
the atomic interaction, the total momentum distribution 
is given as the weighted incoherent superposition of ev¬ 
ery site’s wavefunction’s evolution. Therefore, the whole 
atom cloud could also exhibit a good oscillation behav¬ 
ior with time. However, if the coherence term exists, the 
result would be different, and the collective oscillation is 
destroyed as shown by the simulation in Fig. 3 of main 
body. 


FIG. 6: Evolution of different site’s wave function for 60ms 
with 5Er lattice depth. The simulation includes totally 301 
lattice sites. The wave function (a-f) is initially distributed 
in the —25th, —20th, —15th, —10th, —5th or 0th site respec¬ 
tively. 
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FIG. 7: Evolution of different site’s wave function for 60ms 
with 10Er lattice depth. Other conditions are the same as 
Fig. 3. 
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